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Abstract — The study is devoted to radial pulsations in Population I Cephcids with masses 
from 5.4M 0 to 6M 0 . Solution of the equations of radiation hydrodynamics and turbulent 
convection for nonlinear stellar pulsations was obtained with initial conditions as the core 
helium burning models of computed evolutionary sequences. For each value of the initial mass 
we considered stellar evolution from the zero age main sequence to central helium exhaustion 
with three values of the convective overshoot parameter: a ov = 0.1, 0.15 and 0.2. Models for the 
Cepheid SZ Tau with pulsation period II = 3.149 day can be constructed only for a ov = 0.1 and 
cc 0 v = 0.15. The star is at the evolutinary stage of the second crossing of the instability strip 
and pulsates in the first overtone just near the boundary between the fundamental mode and 
the first overtone. The mass, radius and age of the star are in the ranges 5.46 < M/M 0 < 5.75, 
41.5 < R/Rq < 42.3 and 6.9 x 10 7 yr < t ev < 8.0 x 10' yr, respectively. Predicted period 
change rates are of n — 0.4s/yr. 
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INTRODUCTION 


The pulsating variable star SZ Tau belongs to a small group of Cepheids that are members 
of open clusters. The first report that the Cepheid SZ Tau observed in the extended corona of 
NGC 1647 belongs to this cluster was made by Efremov (1964) on the basis of similar values 
of the distance and the radial velocity of the Cepheid and the stars of the cluster. Later this 
conclusion was confirmed by other observational data (Gieren 1985; Turner 1992; Anderson et 
al. 2013). The study of the Cepheid which belongs to the stellar cluster and determination of 
its fundamental parameters (the mass, radius and luminosity) by methods of stellar pulsation 
theory is of great importance since provides us with an independent method of evaluation of 
the distance and the age of the stellar cluster. 

In the General Catalogue of Variable Stars (Samus et al. 2012) the star SZ Tau with a period 
II = 3.14873 day is referred to Cepheids of type DCEPS with nearly symmetrical light curves 
of small amplitude. Turner (1992) estimated that the average absolute magnitude of SZ Tau 
corresponds to radial pulsations in the Erst overtone. Later this assumption was confirmed by 
Fourier analysis of the radial velocity curve (Antonello and Aikawa 1995). Another confirmation 
of the first overtone pulsation was provided in the work by Bono et al. (2001) where the mean 
stellar radius R = 45.6i? 0 obtained from the radial velocity observations (Bersier et al. 1994) 
was compared with empirical period-radius relation of Calactic Cepheids. It should be noted 
that the mean radius of SZ Tau was evaluated by the Baade-Wesselink method in a number 
of studies (Gieren 1985; Moffett and Barnes 1987; Laney and Stobie 1995; Ripepi et al. 1997; 
Rastorguev and Dambis 2011) but the large spread in values of the radius (34i? 0 < R < 57 R & ) 
is a serious handicap for determining the stellar mass from the period-mean density relation. 

In contrast to most of the Cepheids the light curve of SZ Tau reveals the absence of strictly 
repetitive light variability and this feature substantially complicates the detection of period 
change due to stellar evolution. In particular, Szabados (1977) and Evans et al. (2015) found 
significant erratic variations in the O — C diagram which do not allow for an unambiguous 
detection of the secular period change. However Berdnikov and Pastukhova (1995) showed 
that the O — C diagram can be fitted by the quadratic function and therefore SZ Tau is on 
the evolutionary stage of the second crossing of the Cepheid instability strip and its period 
decreases with a rate of n = —0.49 s/yr. 

The goal of the present study is to theoretically estimate the fundamental parameters of 
SZ Tau and to validate the conclusion of Berdnikov and Pastukhova (1995) on the evolutionary 
status of the star. To this end we carry out the consistent calculations of stellar evolution 
and nonlinear stellar pulsations and determine the conditions for pulsations with period of 
n = 3.149 day, whereas the theoretical estimates of the period change rate are compared with 
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II = —0.49 s/yr. In this work we do not use the period-mean density relation and therefore 
evaluation of the stellar mass does not suffer from uncertainties in the stellar radius estimates 
obtained by the Baade-Wesselink method. Earlier the author employed such an approach for 
determination of the fundamental parameters and the age of the Cepheid a UMi (Fadeyev 
2015). 

Basic equations and methods of their solution are described in previous papers of the author 
(Fadeyev 2013a, 2014, 2015). Results of computations presented below are obtained for the 
initial composition X = 0.7, Z = 0.02, where X and Z are the mass fractional abundances 
of hydrogen and of elements heavier than helium, respectively. Convection in evolving stars 
is treated according to the mixing length theory (Bohm-Vitense 1958) with a mixing length 
to pressure scale height ratio a a = 1.6. Enlargement of the convective core due to convective 
overshooting is taken into account by the parameter a ov = Ar/Hp, where Ar is the radial 
extention of the outer boundary of the convective core and Hp is the pressure scale height. 
The evolutionary computations were done for a ov = 0.1, 0.15, 0.2 and the parameter a ov was 
assumed to be the same from the main sequence to core helium exhaustion. 

RESULTS OF COMPUTATIONS 

To determine the ranges of the initial stellar mass Mzams and the parameter a ow for which 
the radial pulsations with a period of II = 3.149 day become possible during the second crossing 
of the instability strip we computed the grid of evolutioinary tracks for stars with masses 
5.4M 0 < Mzams A 6M 0 . Together with computation of stellar evolution we solved the linear 
adiabatic wave equation (see, e.g., Cox 1983) and evaluated eigenfrequencies of the fundamental 
mode and the first overtone. If the evolutionary model with effective temperature 5400 K < 
T e s A 6700 K has the adiabatic period of one of these oscillation modes close to the value 
II = 3.149 day then from 10 to 15 models of the evolutionary sequence were selected and were 
used as initial conditions for solution of the equations of hydrodynamics. 

From results of hydrodynamic computations we primarily determined the edges of the 
pulsation instability strip, that is the age of the star t ev corresponding to the zero growth rate 
of the kinetic energy. Subsequently for each hydrodynamic model we calculated the periods of 
the fundamental mode and the first overtone using the Fourier transform of the kinetic energy 
of pulsation motions. In the most of considered evolutionary sequences the mode switch from 
the fundamental mode to the first overtone occurs while the star crosses the instability strip. 
In such a case it was assumed that the mode switch is abrupt and the change of the pulsation 
period II within the instability strip is described by two continuous functions. Determination of 
the mode switch boundary is discussed in our previous papers (Fadeyev, 2013b; 2015). Within 
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the evolutionary time interval with continuous change of the period the sequence of the values of 
II obtained from hydrodynamic computations is fitted by the algebraic polynomial of the order 
2 < n < 4. The choice of the polynomial order n obeys the condition that the approximation 
rms error cannot exceed 0.1%. 

For Cephcid hydrodynamic models computed in the present study the most short period 
of the fundamental mode is fl « 4.3 day therefore the following discussion of SZ Tau model 
will be confined to consideration of the first overtone pulsators. The existence of the point of 
the evolutionary track with pulsation period II = 3.149 day depends on both the initial stellar 
mass Mzams and the overshoot parameter a ov . Dependence on the initial mass is illustrated 
in Fig. m where the evolutionary tracks in the Hertzsprung-Russel (HR) diagram are shown 
for the models Mzams — 5.5M 0 and Mzams — 5.8M 0 computed with the overshoot parameter 
a ov = 0.15. For the sake of graphical clarity, we present in the figure only the blue loops 
corresponding to the core helium burning. Evolution along the shown tracks proceeds clockwise. 

As soon as a star with an initial mass Mzams = 5.5M 0 crosses the red edge of the instability 
strip (T e g = 5550 K) it becomes a fundamental mode pulsator. The pulsation switches to the 
Erst overtone at the effective temperature T e g = 5840 K and the pulsation ceases at the blue 
edge of the instability strip (T e g = 6300K). As can be seen in Fig. [lj the model Mzams = 5.5M 0 
with pulsation period n = 3.149 day locates just near the boundary between the fundamental 
mode and the first overtone, 

As the initial stellar mass increases the boundary of the mode switch from the fundamental 
mode to the first overtone moves to higher values of T c g and for stars Mzams > 5.8M 0 becomes 
beyond the blue edge of the instability strip. Therefore, as seen in Fig. [Q the star with initial 
mass Mzams = 5.8M 0 remains the fundamental mode pulsator within the whole instability 
strip (5360 K < T c r < 6250 K). The pulsation periods are confined to the interval 4.45 day < 
n < 6.90 day. 

Conditions for the occurence of pulsations with period n = 3.149 day on the evolutionary 
tracks computed with different values of the overshoot parameter a ov are illustrated in Fig. [2] 
for Mzams — 5.6M 0 . As can be seen, the luminosity of the Cepheid increases with increasing 
a ov (mainly due to the larger mass of the convective core during the main sequence stage), 
whereas the maximum effective temperature of the blue loop decreases. As a result, among 
three evolutionary tracks shown in Fig. [2] only one (a ov = 0.15) can be used for construction 
of the SZ Tau model, the star locating just near the boundary between the fundamental mode 
and the first overtone. The evolutionary track computed with a ov = 0.1 is not appropriate for 
the model of SZ Tau because of too short period (n = 2.81 day) of the first overtone at the 
mode switch boundary. On the other hand, for the evolutionary track computed with a ov = 0.2 
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the blue edge of the instability strip is beyond the turning point of the track so that during the 
stage of increasing effective temperature the Cepheid remains the fundamental mode pulsator. 

It should be noted that for all evolutionary tracks (5.4M 0 < Mzams < 6M 0 ) computed 
with a ov = 0.2 the radial pulsations during the second crossing of the instability strip were 
found to be due to the instability in the fundamental mode. Therefore further discussion will 
be confined to results of hydrodynamical computations done with initial conditions obtained 
with a ov = O.f and a ov = 0.15. 

Conditions of existence of the pulsation with period n = 3.149 day can be represented in 
more detail by the dependence of the period of the first overtone on the evolutionary time. The 
plots for several evolutionary sequences computed with a ov = 0.1 are shown in Fig. [3j Each 
plot ends at the blue edge of the instability strip. For the sake of graphical convenience the 
evolutionary time t ev is set to zero at the mode switch boundary. 

For Mzams > 5.9M 0 the star pulsates in the fundamental mode within the whole instability 
strip. On the other hand, for Mzams < 5.6M 0 the period of the first overtone is smaller than 
II = 3.149 day. Therefore, the model of SZ Tan with the overshoot parameter q ov = 0.1 can 
be constructed for the initial stellar mass Mzams ranging from 5.6M 0 to 5.9M 0 . The principal 
cause of the narrow mass interval is the close position of the Cepheid to the mode switch 
boundary. 

Temporal dependences of the first overtone period in the evolutionary sequences computed 
with overshoot parameter a ov = 0.15 show a similar behavior. The only difference in comparison 
with models computed with a ov = 0.1 is that pulsations with period hi = 3.149 day arise in 
stars with lower masses: Mzams = 5.5M 0 and Mzams = 5.6M 0 . The narrow ranges of initial 
stellar masses allowed us to approximately evaluate the fundamental parameters of SZ Tau 
which are given in the table. 

To compare our theoretical results with observations Fig.[4]shows the plots of hi as a function 
of period II for Cepheid models pulsating in the first overtone. Each curve in Fig. [4] represents 
the evolutionary track in the period-period change rate diagram and evolution proceeds as the 
pulsation period hi decreases. The position of SZ Tau on the diagram and its rms error are 
shown on the basis of data given by Berdnikov and Pastukhova (1995). 

CONCLUSIONS 

The Cepheid SZ Tau is located just near the mode switch from the fundamental mode to 
the first overtone. The proximity of the pulsation period If = 3.149 day to the boundary value 
of the period of first overtone is by an order of magnitude 511/11 ~ 10~ 2 . This is the principal 
cause that stellar evolution and nonlinear stellar pulsation calculations lead to the narrow 
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ranges of the fundamental parameters of the Cepheid. It should also be noted that models of 
SZ Tau obtained for the overshoot parameters a ov = 0.1 and a ov = 0.15 give evidence in favor 
of the small and intermediate convective overshooting during the main sequence and the helium 
core burning stages. 

Predicted period change rates obtained in the present study are in a good agreement with 
the observational estimate (II = —4.9 s/yr) by Berdnikov and Pastukhova (1995) and therefore 
confirm a suggestion that SZ Tau is at the evolutionary stage of the second crossing of the 
instability strip. 

Theoretical estimates of the radius of SZ Tau (41.512© < R < 42.3 Rq) given in the table 
show better agreement with measurements by Ripepi et al. (1997). Authors of this work 
employed the improved version of the Baade-Wesselink method (the CORS method) based on 
the more correct relations between the color index and the surface brightness (Caccin et al. 
1981). For two variants of the CORS method Ripepi et al. (1997) estimated the mean radius 
as R = 41.8-Rq and R = 44.8i?0. 

Hydrodynamic computations of nonlinear pulsations in SZ Tau located just at the mode 
switch are accompanied with great numerical difficulties. This is mainly due to the small decay 
rate of the fundamental mode which leads to slow relaxation of nonregular oscillations after 
attainment of the first overtone limit cycle. Absence of strict regularity in the observed light 
and radial velocity curves of SZ Tau might be due to slow decay of the fundamental mode. 
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Models of Cepheids with period 3.149 day 


C^ov 

^ZAMS /Mq 

M/Mq 

Rj Rq 

L/L e 

T eff , K 

ii, s/yr 

t ev , 10 6 yr 

0.1 

5.7 

5.66 

42.0 

1878 

5870 

-0.42 

72.1 


5.8 

5.75 

42.3 

2008 

5950 

-0.43 

69.2 

0.15 

5.5 

5.46 

41.5 

1819 

5860 

-0.38 

80.2 


5.6 

5.56 

41.8 

1968 

5960 

-0.43 

77.0 
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Figure captions 


Fig. 1. Evolutionary tracks of core helium burning stars with initial masses 5.5M 0 and 5.8M 0 
computed with the overshoot parameter a ov = 0.15. Parts of the track corresponding to 
the pulsational instability during the second crossing of the instability strip are shown by 
the dashed and dotted lines for the fundamental mode and the first overtone, respectively. 
The Elled circle indicates the position of the Cepheid pulsating in the first overtone with 
period II = 3.149 day. 

Fig. 2. Same as Fig. [Q but for the initial stellar mass M Z ams — 5.6M 0 with a ov = 0.1, 0.15 and 

0 . 2 . 

Fig. 3. The first overtone pulsation period IT as a function of the evolutionary time t ev during the 
second crossing of the instability strip. Solid lines show evolutionary sequences computed 
for a ov = 0.1. The evolutionary time t ev = 0 corresponds to the mode switch from the 
fundamental mode to the first overtone. 

Fig. 4. The period change rate II as a function of period II for Cepheids pulsating in the first 
overtone. The solid and dashed lines show the evolutionary sequences computed for the 
overshoot parameter a ov = 0.1 and a ov = 0.15, respectively. The values of M z ams are 
given at the curves. The position of the Cepheid SZ Tau corresponds to the observational 
estimates by Berdnikov and Pastukhova (1995). 
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Figure 1: Evolutionary tracks of core helium burning stars with initial masses 5.5M 0 and 
5.8M 0 computed with the overshoot parameter a ov = 0.15. Parts of the track corresponding 
to the pnlsational instability during the second crossing of the instability strip are shown by 
the dashed and dotted lines for the fundamental mode and the first overtone, respectively. The 
filled circle indicates the position of the Cepheid pulsating in the first overtone with period 
IT = 3.149 day. 
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log L/L 



Figure 2: Same as Fig. [T]but for the initial stellar mass M Z \ms — 5.6M 0 with ct ov = 0.1, 0.15 
and 0.2. 
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Figure 3: The first overtone pulsation period II as a function of the evolutionary time t ev 
during the second crossing of the instability strip. Solid lines show evolutionary sequences 
computed for a ov = 0.1. The evolutionary time t ev = 0 corresponds to the mode switch from 
the fundamental mode to the first overtone. 
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Figure 4: The period change rate 14 as a function of period n for Cepheids pulsating in the 
first overtone. The solid and dashed lines show the evolutionary sequences computed for the 
overshoot parameter a ov = 0.1 and a ov = 0.15, respectively. The values of Mzams are given at 
the curves. The position of the Cepheid SZ Tau corresponds to the observational estimates by 
Berdnikov and Pastukhova (1995). 
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